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Abstract 

The charge and baryon balance functions are studied in the co- 
alescence hadronization mechanism of quark-gluon plasma. Assum- 
ing that in the plasma phase the qq pairs form uncorrelated clusters 
whose decay is also uncorrelated, one can understand the observed 
small width of the charge balance function in the Gaussian approx- 
imation. The coalescence model predicts even smaller width of the 
baryon-antibaryon balance function: a^g/aj^ = ^2/3. 



1. Introduction 



It was recently argued jT] that the constituent quark coalescence model j2] 
can explain the small pseudo-rapidity width of the charge balance function 
O E] observed in central collisions of heavy ions at RHIC , as compared 
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to that measured in pp collisions The charge (or baryon number) bal- 
ance function measures the relative distance between produced positive and 
negative charges. Encouraged by the success of this simple argument, in the 
present paper we apply the same idea to the balance function measuring the 
baryon-antibaryon compensation. 

Our aim is to discuss correlations between produced hadrons, and thus 
the quark coalescence model must be generalized to consider possible pre- 
existent correlations inside the system. In a state rich in quark- ant iquark 
pairs arising from initial (charge and flavor neutral) gluon fragmentation we 
can supplement the standard formulation of the model by the hypothesis 
that -before hadronization- the QGP forms charge- and flavor-neutral qq- 
clusters, which seems natural in a gluon-dominated plasma. In the final, 
pre-hadronization stage the clusters dissociate, quarks and antiquarks coa- 
lesce into observed hadrons while gluons required to balance color charge 
form additional neutral quark clusters (and again gluons) and the process is 
repeated till hadronization is completed. 

For this reason we can assume for the purpose of evaluating the balance 
correlation that at the last stage before hadronization the plasma consists 
of isotropic qq clusters which are charge-, flavor- and (though natural, here 
important) baryon number- neutral. In the coalescence model a charged 
meson is formed from two such qq clusters. The charged bosons final state 
balance function is thus reduced by a factor 1/ V2 with respect to the width 
of source distributions, since as is well known, the width of the distribution of 
the average of any two independent variables (g, q in different qq) is reduced 
by the factor 1 / ^/2 with respect to the width of their own distributions. This 
was the argument in Ref. pp. 

Here we observe that since baryons are formed from three constituent 
quarks, we expect in this case a reduction by a factor l/y/S, and thus an 
even further reduction (by the factor ^2/3 beyond charge balance narrow- 
ing) of the width of the baryon balance function. A particularly interesting 
case is the balance function of neutral strange baryons {A{sud)) since the co- 
alescence formation requires here under any possible microscopic mechanism 
the participation of three independent {ss, uu, dd) neutral quark clusters. 

The confirmation of this prediction by the future data could be consid- 
ered as an important indication that the coalescence model is the dominant 
mechanism in the hadronization of the quark-gluon plasma. The coalescence 
model is the key ingredient required to describe the relative yields of hadrons 
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within the statistical hadronization approach jB]: the probabihty of produc- 
tion of a hadron is by assumption a product of probabihties of finding the 
constituent quarks of the hadron considered. 

With this formulation all consequences of the coalescence model concern- 
ing single particle distributions are of course unchanged. In particular, the 
spectacular successes of the model in description of single particle spectra 
[21 n\ and flow parameters [S] in the central collisions of heavy ions remain 
valid. Thus our extension is a natural attempt to make the coalescence model 
more complete. 

To proceed, we need to define some properties of the neutral clusters. To 
keep the discussion as simple as possible, we shall assume that (i) a cluster 
decays isotropically in its rest frame, and (ii) momenta of the decay products 
are uncorrelated. These two properties are sufficient to derive most of our 
conclusions. To keep the problem as simple as possible, we shall also assume 
that the clusters themselves are uncorrelated, as would be expected if these 
arise from gluons. 

Using this scenario, we give predictions for the (pseudo) rapidity width 
of the charge and baryon balance functions, including corrections due to the 
finite acceptance region. In the next section „ we introduce the notation. The 
analysis of the charge balance function given in jl] is shortly summarized in 
section 3. Estimate of the baryon balance function is given in section 4. Our 
conclusions are listed in the last section. 

2. Balance functions 

In this paper we shall only discuss the neutral systems, u.e. the systems with 
vanishing total charge and baryon number. In this case the balance function 
of the "charges" denoted by + and — is expressed in terms of the single- 
and double- D particle functions jH 1^ as follows: 



5(A2|Ai) 



D(+,A2|-,Ai)-D(+,A2|+,Ai) 



+ 




+ 



D(-,A2|+,Ai)-D(-,A2|-,Ai) 



./9 



(1) 




Here 




dyidy2 



(2) 



3 



where dui/dy and dnfj/dyidy2 are the corresponding particle densities. 

The measurement of the electric charge balance function performed by the 
STAR collaboration requires both particles to be in a finite rapidity accep- 
tance interval —A < yi,y2 < A while the difference S of (pseudo)rapidities 
is kept fixed. This suggests a change of variables: 

2/1+2/2 _ 

2/1 - 2/2 = d; — - — = z.hx (4) 

With this notation, the integration in Eq. (??) must be restricted to 5 
being fixed, while z must be kept inside the interval —A < z < A where 

A = A - ^-^.lixa (5) 

We shall use this prescription also for the baryon number which is of interest 
in this paper. 

Thus the balance function we discuss is defined as 



/_ A [d'^n^_ I dzdS — d'^n^^ / dzdS] 
i_A dy+dn+/dy+ 

/Ja d^ [(i^^-H- / dzdS — d^n / dzd5] 



lib (6) 



J^^ dy^dn^/dy_ 
where ± refers here to the baryon number. 

3. Balance function in the coalescence model 

We first consider the production of charged pions IJ. In formation of such a 
pair, the dominant process is that involving one ^/-cluster and one D-cluster: 

{uu) + (dd) (ud) + {ud).lml (7) 

The corresponding two particle momentum distribution is: 

p{p+,p-) = J dPidP2pc{Pi)pc{P2) j dpu^dpu^dpd^dpd^ 

f{Pl,Pu,)fiPl,Pu,)fiP2,Pd,)f{P2,Pd,) 

6[p+ - {pu, +pdj/2](5[p_ - {pu, +Pd2)/2] 

Gm{Pm - Pd2)GmiPdi - Pu2)i'm'^ (8) 
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where Pc(-P) is the distribution of clusters, f{P,p) the cluster dissociation 
probability, that is the probability of finding a quark of momentum p in the 
cluster of momentum P, and G{p — p) is the probability of coalescence, i.e. 
that a quark-antiquark pair forms a hadron. The momenta of quarks are 
denoted by p, those of antiquarks by p. 

Since one pair of clusters cannot produce two mesons of the same charge, 
only p{p+,p^) given by (??) contributes to the balance function. In the 
following we shall only consider rapidity and pseudo-rapidity distributions. 
To simplify the technicalities we take all functions in (??) in the Gaussian 
form: 



Pc{Y) = -^e-^V^^ ; f{Y, y) = J-e-(^-)Vc^ ; Im^ (9) 



where is the total number of clusters. It turns out that the balance function 
does not depend on the shape of G{u). 

Introducing (??) into (??) allows to perform all integrations. The result, 
substituted in (??), gives 



One sees that in the limit of large acceptance (A — > oo) the width of the 
balance function is determined by the same parameter c which describes the 
cluster dissociation width. The effect of finite acceptance. A, is seen as the 
ratio of the two error functions. Note that in the Gaussian distribution model 
the balance function is entirely independent of the coalescence probability 
distribution G{u). 

4. Baryon balance function 

The dominant process in the production of baryon- antibaryon pairs involves 
three clusters: 




(10) 



{qq) + {qq) + {qq) 



(qqq) + {qqq)lbl 



(11) 



Consequently, the distribution of the BB pair is: 




dPidP2dP3Pc{Pi)pc{P2)Pc{P3) / dpidp2dp2,dpidp2dp2, 
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f{Pl,Pl)f{P2,P2)f{P3,P3)f{Pl,Pl)f{P2,p2)f{P3,p3) 
S[p - {pi +P2+ P3)/3]5[p - ipl +p2+ P3)/3] 

GsiPl -P2,P2 -P3,P3-Pi)Gb{Pi -p2,p2 -p3,p3-pl)- W 

Introducing Eq. (??) into Eq. (??), performing the integrations and putting 
the result into Eq. (??), one obtains: 

One sees that in the hmit of large acceptance the width of baryon balance 



function is again determined by the parameter c but it is by the factor y2/3 
smaller than that for obtained for the charged pion balance functions given 
by Eq. (??). The corrections due to finite acceptance are also given by ratio of 



two error functions but their arguments are multiphed by y3/2, as compared 
to those in Eq. (??). Therefore they are slightly less important than in the 
pion case. 

It is natural to expect that the {qq) cluster dissociation probability is 
isotropic. As discussed in JT] , this puts a strong constraint on the parameter 
c responsible for the cluster decay distribution. Indeed, the isotropic decay 
implies that the decay distribution in pseudo-rapidity is 

corresponding to the (pseudo) rapidity distribution of statistically produced 
(massless) particles. The width of this distribution is (|?7|) = / d7]\7]\p{7]) = 
log 2 ^ 0.69. In order to approximate the distribution (??) by a Gaussian of 
the same width we must use: 

c = V^log2 ^ 1.23./nl (15) 

Although Eq. (??) refers to pseudo-rapidity distribution, we shall take 
this value of c to describe as well the rapidity distribution in cluster decay^. 
As shown in jT], this value is compatible with the STAR data [H], provided 
that the finite acceptance corrections are included {A = 3.5 was assumed). 

"'^One can show that the parameter c of the rapidity distribution cannot be larger that 
that describing the pseudo-rapidity distribution. 
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These corrections are important: they reduce the width from 0.69 to 0.55, 
compatible with data. 

The comparison of the resuhs for charge and baryon balance functions 
is shown in Fig. ^ where the expected widths, calculated from Eq. (??) and 
Eq. (??) with c given by Eq. (??) and A = 1.3, are plotted as function of 
the cluster distribution width A, to show (in)sensitivity with regard to the 
cluster distribution. Also in Fig. Qthe ratio of the charge (m: meson) and 
baryon balance function widths is shown. 
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Figure 1: Widths of balance functions versus the width A of the pre- 
hadronization cluster rapidity distribution. 

As discussed in Ref. P, {\y\) for pions agrees reasonably well with the 
STAR data in a rather broad range of A. One also sees that the ratio of 
baryon and pion widths is predicted to be .87 and is insensitive to A. This 
value is larger than y^2/3) ~ .82, expected for large rapidity acceptance. 

In Fig. 12 the same quantities are plotted versus A, the size of the accep- 
tance region for A = 3.5. One sees that the widths change rather rapidly 
when A varies from 1 up to 3. Starting from A = 4 the effect of acceptance 
practically disappears. 



5. Comments and conclusions 

Several comments are in order. 
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Figure 2: Widths of balance functions versus the acceptance A. Vertical 
dashed line indicates the STAR acceptance. 

(i) We have ignored the transverse motion of the clusters (the distribu- 
tion (??) is only valid for clusters with vanishing transverse momentum). 
The effect of transverse motion is expected to reduce the decay width in ra- 
pidity ^U]. This effect was discussed in Ref. P] and shown to be small for 
transverse velocities up to 0.4 . It may be interesting, however, to investi- 
gate -both theoretically and experimentally- the balance functions at larger 
transverse momenta. Such studies may provide additional information about 
the properties of the transverse flow. 

(ii) We have also ignored the difference between rapidity and pseudo- 
rapidity. Therefore, our predictions of the balance function width in rapidity 
could be somewhat overestimated. The width ratio between baryons and 
pions is, however, much less sensitive to this effect and thus remains a much 
more solid prediction. 

(iii) The correlations in cluster decay were neglected. As shown in P, 
introducing correlations may change somewhat the width of the balance func- 
tions. This may be an interesting point for future research. 

(iv) As we were using Gaussian approximation, we cannot discuss se- 
riously the shape of the balance function. Although the present data are 
consistent with Gaussian shapes, it would be interesting to investigate the 
problem more thoroughly when better data become available. 

(v) We have only considered the dominant processes leading to meson 
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and baryon formation. The corrections from sub-dominant ones give, gener- 
ally, contributions to balance functions which are much broader than those 
described in the present paper. 

(vi) The effect of resonance decay was entirely neglected. This serious 
problem seems, however, beyond reach of our approach. In this context, it 
should be emphasized that the argument presented in this paper applies to 
any kind of baryon-antibaryon pairs. Thus we feel that measurements of 
(multi) strange baryons may be particularly interesting, because -due to their 
large mass- they are less sensitive to corrections related to resonance decays 
and transverse motion. 

In conclusion, we have calculated the width of the charge and baryon 
number balance functions in the coalescence model generalized to include 
possibility of correlations between constituent quarks and antiquarks in the 
form of charge- and flavor- neutral, isotropically dissociating clusters. The 
existing data on charge balance functions are explained. The width of the 
baryon balance function at large acceptance is predicted to be smaller by 
nearly a factor ^2/3. Finite acceptance corrections were estimated. We be- 
lieve that future measurements of the baryon balance functions may provide 
an interesting test of the quark coalescence hadronization model. 
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